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ABSTRACT 


A method  is  developed  for  determining  the  proportion  of 
parcels  which  "overflow"  in  an  automatic  sorter^  in  terms  of 
the  capacities  of  the  buffers  employed,  the  number  of  operator 
lanes,  and  other  system  parameters » The  treatment  extends 
preA^ious  ones  by  permitting  consideration  of  arbitrary  buffer 
sizes,  and  by  yielding  time-varying  as  well  as  steady  state 
solutions..  Numerical  results  for  the  steady  state  case  are 
given,  and  display  a drastic  increase  in  overflow  rate  when 
the  number  of  operator  lanes  exceeds  a critical  value. 


ANALYSIS  OF  OVERFLOW  RATE 


IN 

A SORTING  SYSTEM 

by 

B.K.  Bender  and  A«J.  Goldman 
INTRODUCTION 

This  report  treats  a question  arising  in  connection  with  the 
parcel  post  sorter  being  designed  for  the  Post  Office  Department 
by  Food  Machinery  and  Chemical  Corporation  (FMC)»  A queueing 
approach  to  the  problem^  assiiming  random  (Poisson)  arrival  and  service 
times_,  was  discussed  in  Volume  9 of  FMC's  3/31/6I  Gateway  Project  Report 
Under  the  assumptions  of  uniform  arrival  and  service  rates,  plus  a 
hj^othesis  of  equal  distribution  of  parcels  among  cross  belts,  the 
"steady- state  no-buffer"  case  was  solved  by  John  Weaver  of  FMC  in 
the  9/26/60  Technical  Memorandum  No»  STM-8,  "Analysis  of  the  Jamming 
Phenomenon  in  the  Matrix  Sorting  Machine". 

Using  the  same  "iiniformity"  assumptions  as  in  the  Technical 
Memorandum  cited,  we  shall  describe  a solution  method  valid  for 

(a)  arbitrary  buffer  sizes, 

(b)  the  t Ime- varying  case  as  well  as  the  steady- state  case,  and 

(c)  arbitrary  distribution  of  parcels  among  cross  belts. 

The  generality  of  the  solution  is  perhaps  greater  than  seems  necessary 
now,  but  may  prove  useful  in  future  work.  It  seems  appropriate  to 
point  out  explicitly  that  the  methods  employed  are  purely  anal;>’t;ic 
(i.e.,  no  simulation  is  required)  and  also  mainly  elementary;  only 
simple  algebra  and  probability  theory,  together  with  mathematical 
induction,  are  used. 
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A brief  description  of  the  problem  is  given  in  Section  2^ 


while  Section  3 contains  a detailed  account  of  the  "mathematical  model" 
to  be  employed.  The  presentation  of  the  general  solution  method  is 
begun  in  Section  h,  but  the  discussion  is  then  narrowed  to  a special 
"steady  state"  situation  believed  to  be  of  particular  interest.  The 
resulting  formulas,  which  represent  a generalization  to  arbitrary 
buffer  size  of  those  given  in  FMC's  memo,  are  used  to  derive  numerical 
results  for  realistic  values  of  the  of  the  relevant  parameters.  These 
results  show  that  the  overflow  rate  increases  drastically  when  the  number 
of  operator  lanes  reaches  a certain  level  (determined  by  the  parameters), 
so  that  additional  lanes  would  accomplish  relatively  little  sorting.  The 
reader's  attention  is  especially  directed  to  this  material.  Numerical 
results  for  a wider  range  of  cases  are  being  calculated,  and  will  be 
reported  separately. 

In  Section  5 the  description  of  the  general  solution  method  is 
resumed  and  completed.  The  method  is  illustrated  by  three  applications. 
The  first  sketches  the  details  for  treating  a small  buffer,  with  room 
for  only  one  parcel.  The  second  refers  to  a steady  state  situation  in 
which  cross  belts  move  faster  than  operator  lanes.  The  third,  given 
in  Section  6,  extends  FMC's  "no  buffer"  analysis  to  the  case  in  which 
the  distribution  by  destinations  of  parcels  varies  with  time  during 
the  sorting  system' s operation. 
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2.  STATEMENT  OF  THE  PROBLEM 


Concept-ually^  each  parcel  arriving  at  the  sorter  is  "coded" 
by  an  operator  and  placed  on  a moving  belt  called  an  operator  lane» 

This  belt  passes  over  a series  of  cross  belts^  also  moving*  At  the  end 
of  each  cross  belt  is  a set  of  bins  for  accumulating  parcels  to  a 
partlcnlar  set  of  destinations.  When  the  parcel  moving  along  the 
operator  lane  passes  over  the  proper  cross  belt;,  it  is  transferred  to 
that  cross  belt  and  eventually  to  the  proper  destination  bin. 

In  practice^  however^  parcels  are  coded  not  by  one  but  by  any 
of  several  operators^  and  are  placed  on  any  of  several  operator  lanes^ 

.all  passing  over  the  same  set  of  cross  belts.  The  complication  arises 
■when  a parcel  Intended  for  a particular  cross  belt  reaches  that  belt 
and  finds  the  passing  section  of  it  to  be  already  occupied  by  another 
parce.l  (from  an  "earlier"  operator  lane).  The  new  parcel  cannot  occupy 
that  same  section  (assuming  this  were  physically  possible)^  since  the 
code  ca,nn,ot  be  retained  for  two  different  destinations;  one  parce,l  woul.d 
end  in  the  wrong  bln. 

One  possibility  is  to  ha.ve  each  parcel^  for  which  space  on  the 
appropriate  cross  belt  is  not  available  when  needed.,  simply  continue 
to  the  end  of  the  operator  .lane  and  then  be  rehandled  or  recycled.  This 
is  the  "no  buffer"  case  mentioned  above.  An  alternative  is  to  place 
buffers  at  the  Intersection  points  of  operator  lanes  and  cross  belts^ 
so  that  a parcel,  could  wait  in  the  appropriate  buffer  until  space  for  it 
became  available  on  the  passing  cross  belt.  Since  definite  costs  are 
associated  with  both  the  addition  of  buffers  a.nd  the  recycling  of  parcel-S;, 
it  was  desired  to  ascertain  by  how  much  Installation  of  buffers  of  vario’as 
capacities  would  reduce  the  fraction  of  parcels  requiring  3:*ecycling_,  and 
how  this  fraction  depends  on  the  number  of  operator  lanes. 

A somewhat  different  technique  may  occur  to  the  reader.  One  would 
place  bins  for  an  especially  frequent  destination  at  the  ends  of  several 
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cross  "beltsj  if  a parcel  for  such  a destination  was  unable  to  transfer 
to  the  first  appropriate  cross  belt,  it^  would  have  "another  chance" 
further  along  the  operator  lane.  This  would  involve  certain  equipment 
complications;  what  is  done  instead  is  to  Install  bins  for  major 
destinations  directly  "on"  the  operator  lanes,  between  the  operators^ 
coding  stations  and  the  system  under  analysis  here* 
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3.  THE  MATHEMATICAL  MODEL 


Our  mathematical  model  of  the  situation  is  a discrete  one.  A 
belt  section  is  in  effect  identified  -with  some  appropriate  point  on 
that  section;  the  parcel  (if  any)  carried  by  the  section  is  thought  of 
as  concentrated  at  the  corresponding  pointy  and  the  section  is  considered 
to  pass  another  belt  at  the  mopnent- the  corresponding  point  passes  it. 

f 

This  idealization  simplifies  the  description  of  the  analysis,  but  does 
not  appear  to  restrict  its  validity.  It  also  permits  us  to  think  of 
time  as  moving  in  discrete  steps. 


Some  more  specific  assumptions  vill  nov  be  listed.  The  second  of 
them  could  apparently  be  considerably  weakened  (e.g.,  k^^  instead  of  k^') 
without  causing  much  extra  conceptual  difficulty,  thou-gh  the  details  of 
the  derivations  wo'uld  become  substantially  more  complicated.  The  notation 
(i,j)  will  be  used  for  the  intersection  of  the  i-th  operator  lane  and 
the  j-th  cross  belt. 


(l)  Uniform  Motion.  Sections  of  the  i-th  operator  lane  arrive  at 
(i,j)  at  a constant  rate,  and  the  same  holds  for  the  j-th  cross  belt. 
Tlius  we  can  define 


t^j=time  between  arrivals  of  successive  sections  of  operator  lane 
at  (i,j), 

( g) 

t. /=tlme  between  arrivals  of  successive  sections  of  cross  belt 
ij 

at  (i,j)o 


(2)  Synchronlzatldno  This  assumption  has  three  parts.  First,  for 
each  (i,j)  with  i > 1 there  is  a number  c. . > 1 such  that 

ij  - 

(c)  ! 

c.  ,t.  / =time  between  the  arrival  of  a section  of  cross  belt 
ij  ij  I 

at  (i-l,j)  and  the  same  section's  arriA’^al  at  (i,j). 

For  constant  belt  speed,  this  says  that  (l-l,j)  and  (i,j)  are  c^^ 
"section  lengths"  apart  along  the  j-th  cross  belt*  Second,  there  is  an 
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integer 
lane  speed, 


representing  the  ratio  Of  cross  belt  speed  to  operator 
such  that 


t . =k,t^?^ 
IJ  J Ij 


Third,  at  some  initial  time  it=t*.,  ah’ operator  lane  section  and  a 

cross,  belt  section  arrive  simultaneously  at  (l,j).  From  this  and 

the  second  assumption  it  follows  that  every  arrival  of  an  operator 

lane  section  at  (l,j)  is  simultaneous  with  the  arrival  of  a cro§s 

belt  lane  at  (l,j),  but  not  conversely  if  k.  >1* 

J 


One  consequence  of  assumption  (2)  is  that  all  times  of  interest 
in  connection  with  (i,j)  will  be  of  the  form 

f c ) 

t=t*.  + Nt^.{  , N=0,l,2, 

a cross  belt  section  will  pass  (i,j)  at  each  of  these  times,  but  an 

operator  lane  section  will  pass  if  and  only  if  N is  a multiple  of  k.« 

J 

To  simplify  notation,  we  adopt  the  abbreviation 


(3) 


, (^)  (c) 

trV=  t*.  + wt)/ 

ij  ij 


Independence « Suppose  we 


set 


p^j(t)=probabillty  that  at  time  t a parcel  destined  for  tha 
j-th  cross  belt  arrives  at  (i,j)  along  the  1-th 
operator  lane. 


We  assume  that  the  event  defining  p. .(t)  is  independent  of  both  (a) 

1 J 

the  contents  of  the  buffer  at  (l,j)  at  time  t,  and  (b)  the  state 

(occupied  or  unoccupied)  of  the  section  of  the  j-th  cross  belt  (if  any) 

arriving  at  (i,j)  at  time  t.  In  effect,  we  assiime  that  the  destination 

of  a particular  parcel  entering  the  sorter  is  Independent  of  that  of 

any  other  parcel,  and  also  of  the  choice  of  an  operator  lane  for  the 

parcel.  This  topic  is  discussed  in  the  FMC  Technical  Memorandum  qited 

earlier  ("Step  V"),  Note  that  p.  .(t[^^)=0  unless  N is  divisible-=by  k., 

IJ  J 
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The  mathematical  model  will  now  he  described  in  more  detail^  by 
specifying  what  may  happen  at  any  intersection  (i^j)  at  any  of  the 
(discrete)  times  t concerning  us.  These  times  are  those  at  which  sections 
of  cross  belt  arrive  at  (i.j).  i*e.  t^^^,  = tf.  + N for  N=0,1.2. etc. 

Suppose  first  that  is  not  a multiple  of  so  that  no  sfsction 
of  the  i-th  operator  lane  also  arrives  at  time  t.  The  arriving  section 
of  the  j-th  cross  belt  may  either  be  full  (i»eo,  already  occupied  by 
a particle)  or  emptyj  the  latter  is  always  the  case  if  i=l  (l.e.^  for 
the  first  operator  lane).  If  the  buffer  at  (l;,j)  is  empty^  or  the  passing 
section  of  the  cross  belt  is  full,  then  no  change  ensues.  But  if  the 
buffer  is  not  empty  and  the  passing  section  ^ empty  then  one  of  the 
parcels  in  the  buffer  (it  doesn't  matter  which,  for  oiur  purposes)  is 
transferred  to  the  passing  section.  In  this  connection  it  is  convenient 
to  define 


E. .(t)  = probability  that  the  section  of  cross  belt  passing 

X J 

(i,j)  at  time  t is  empty. 


Next  suppose  that  N ^ a multiple  of  k . , If  the  arriving  section 
of  the  operator  lane  carries  no  parcel  for  the  j~th  cross  belt,  and 
the  passing  section  of  the  choss  belt  is  full,  then  no  change  ensues. 
The  same  is  true  if  no  parcel  for  the  cross  belt  arrives,  the  passing 
section  of  the  cross  belt  is  empty,  and  the  buffer  at  (i, j)  is  empty. 
If  the  buffer  is  non-empty  and  a space  on  the  cross  belt  is  available 
then  one  oJ^  the  parcels  in  the  buffer  is  transferred  to  the  cross  belt; 
the  new  parcel  (if  one  arrives)  goes  into  the  buffer.  If  a new  parcel 
arrives  to  find  an  empty  buffer  and  an  empty  passing  section  of  cross 
belt,  it  is  transferred  directly  to  the  cross  belt.  To  describe  the 
remaining  cases,  it  is  convenient  to  define 


B 


b(l,j)  = capacity  of  the  buffer  at  (i,j), 

. (t)  = probability  that  buffer  at  (i,j)  contains  exactly 

1 jx 

k parcels  at  time  t. 
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so  that  B.  = 0 for  k > Assume  that  a new  parcel  for  the 

j-th  cross  belt  arrives,  but  that  the  passing  section  of  this  belt  Is 
occupied.  If  the  buffer  at  (i,j)  contains  fewer  than  b(i,j)  parcels, 
the  new  parcel  simply  goes  into  the  buffer.  But  if  the  buffer  is  full 
(i.e.,  contains  b(l,j)  parcels),  then  the  new  parcel  counts  as  overflow 
and  disappears  so  far  as  further  analysis  Is  concerned.  We  set 

V.  .(t)  = probability  that  an  overflow  at  (l,j)  occurs 
at  time  t. 

Finally,  to  simplify  notation  we  shall  generally  abbreviate 


if.  SOLUTION  OF  THE  PRINCIPAL  CASE 


Since  the  performance  of  the  sorter,  in  the  present  context, 
is  measured  hy  the  rate  of  overflow,  determination  of  the  functions 
. (t ) is  our  main  objective.  The  preceding  description  of 
the  mathematical  model,  in  particular  the  Independence  Ass^umptlon, 
enables  us  at  once  to  write  the  equation 


(1) 


IJ  IJ 


[l-E, 


(N) 

ij 


] * 


Regarding  the  functions  p(^^  as  "known",  we  see  that  eq(l)  penults 

( n) 

the  calculation  of  the  functions  V. . if  we  can  somehow  determine 

f n''  ( n) 

the  functions  E.  / and  B.  especially  for  k=b(i,j).  A method 

for  determining  these  two  families  of  functions  will  now  be  described. 


The  method  Involves  a fixed  value  of  j,  i.e.  a single  cross-belt 
and  proceeds  by  Induction  on  the  operator  lane  index  1.  For  i=l 
(referring  to  the  operator  lane  met  "first"  by  every  cross  belt)  we 
clearly  have 


(2)  bW=1, 


Bt^,^  = 0 for  k > o 
Ijk 


(3) 


= 1. 


Of  course  no  buffer  would  be  needed  at  (l,j),  i.e.  b(l,j)=o. 


For  the  "induction  step",  suppose  that  the  functions  E; 


(N) 


i-1,  j 


M 


and  B'U'.(  have  been  determined:  we  must  show  how  to  find  the  functions 
1-1,  jk 

and  b|?^.  Since  the  section  of  the  j-th  cross  belt  which  passes 
(i,j)  at  time  is  the  same  one  that  passed  (l-l,j)  at  time 

we  have 

(4)  e(“)=E.  , .(t!’?l  = ,.t(  = ))4..(tW-c.,t(!lB.  , . 

IJ  1-1^ IJ  IJ  IJ  ij  IJ  IJ  1-1, jo'  Ij  IJ  ij 
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This  shows  how  to  determine  the  function  from  the  known  function^ 

0^^)  = and 

(¥e  note  in  passing  that  (2)^  (3)  and  (4)  admit  the  somewhat  more 
explicit  solution 


(5)  E?"!  j,  £ e B . (t«-  | 

-{La  s=r  r=2  s=r+l 


c 


for  i > 1^  where  the  second  product  is  taken  as  unity  for  i=2;  this 

follows  by  Induction  on  i«)  The  remaining  task  is  to  show  how 

the  functions  can  be  determined  from  the  functions  E^^^  and  p(^l 

ijk  ij 

For  this  purpose^  we  examine  what  possible  numbers  of  parcels  in 
buffer  (i^j)  at  time  t^^^,  and  what  intervening  events^  could  lead 

(n+i) 

to  the  bujffer  containing  exactly  k parcels  at  time  t.  , o Using 

^ J 

the  detailed  description  of  the  model  given  in  the  last  section, 
we  obtain  the  following  equations: 


,(N)  _ M r (N)_^  MMi  ^ M p(N)  (N) 


for  k=0. 


(6)  B^^y  = B [qi.^+  Eyypyy]  + B^y  Eyyqy: 

IJO  IJO  IJ  ^IJ  Ijl  IJ  ^IJ 

(7)  b!":")  . b!*?)  , . B?")  [1")  (I-e!"'))  + 

Ijk  ij,k-l  -^IJ  IJ  ijk  ^j  ^ IJ  ^ Ij  IJ 

for  0 < k < b(i,j), 


B b!™)  ,,  E^”) 

ij,k+l  ^ij  ij 


(8)  Biri). ,,  = b'?, 


(N)  j.(K)] 


+ p: . 

ij  ij 


ij,b(i, j) 
for  k=b ( 1 , j ) . 


ij 


In  addition,  the  obvious  Identity 

B{1, j) 


(9)  £ 

k=o 


B?«  = 1 

ijk 


expresses  the  fact  that  the  number  of  parcels  in  the  buffer*  at  (l,j) 
is  always  some  Integer  between  0 and  the  buffer's  capacity,  inclusive 
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So  far  the  analysis  has  been  carried  out  on  a rather  high  level 
of  generality*  The  remainder  of  this  general  analysis  will  be  deferred 
to  Section  5^  and.  attention  will  now  be  focussed  on  a more  special 
situation  of  particular  Interest^  which  we  might  call  the  principal  case* 
The  first  hypothesis  delimiting  this  case  Is 
(10)  = 1, 

l*e«  equal  speeds  for  operator  lanes  and  the  J-th  cross  belt*  The  second_, 
Involving  a kind  of  "homogeneity"^  Is  that 

(lOa)  = Cj  = integer,  tjj)  = t*.  = t|, 


a'll  Independent  of  1«  The  third  assumption  Is  that^  to  a satisfactory 
approximation;  the  steady  state  Input  condition 

(,ll)  p,  .(t^^^)  = p.^^^  = Independent  of  N;  for  all  1 

ij  ij 

holdsj  this  might  occujr  If  the  fluctuations  In  p^^^  are  small  enough  to  be 


Ignored  over  the  lorg  run*  The  same  Is  then  true  of 


- 1 


J , O I 


Under  these  hypotheses;  It  Is  known  (a  few  more  details  are  given  at 
the  end  of  Section  5)  that  the  j-th  cross  belt  will  grad.ually  settle 
down  so  that  the  functlo.ns  and  for  all.  operator  lanes  "1"; 
obey  conditions 


(12)  = E 


(n) 

. . = B_.  ; independent  of  N 


analogous  to  ( 1.1 ) 


/ jq-  \ 

The  determination  of  the  functiohs  B:  (which  are  now  constants  B,  =,  ) 

ijk  ^ ijk- 

can  be  carried  out  especially  simply  for  -^e  principal  case*  Since 

= B^^,^  and  = E^^%  eqs  (6)  and  (j)  can  be  rearranged  to  yield 

Ijx  IJK  ij  ij 

(13)  B . „ q.  „ E,  . = B . . p.  . [l-E, , ] . 

Ijl  ^j  ij  ijo  ^ij  Ij”"  ^ 


®lj;k+l  ‘^.j  ®lj  ®ljk  \j  ^ “ \l,k~l 

for  o < k < b (l; j)o 
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From  this  we  can  prov®,  hy  Irtduetlon  on  that 

(^5)  9ij  Eij  = Pij  [l-Eij]  for  o < k < b(l,j). 

For  k=o,  eq  (I5)  is  precisely  (13)*  If  eq  (15)  holds  vith  k replaced 
by  k-l_,  then  eq  (l4)  yields 


®ij,k+l  ®ljk  ^ “ ®ijk 


Ij 


" ®i  jk  j j 


so  that  (15)  also  holds  for  k*  This  completes  the  induction  proof, 
If  now  define  an  auxiliary  quantity 


(16)  t.. 


ij 


i.J 


then  it  readily  follows  f^ai  #q  (15)  that 


hjk  = (hj)  hjo 


for  o<k<b  (l^j). 
Substitution  from  (17)  into  (9)  yields 


1 = E 

k=o 


k=o 


ij'  ijo  ijo  - ' ij' 


so  that 

(18)  B,  . = 1/  s 


ijo 


k=o 


ij' 


ij 


where  the  first  formula  must  be  used  If  = 1^  and  is  preferable  for 
numerical  work  if  v.  . is  near  unity*  As  a check._,  the  combination  (17) 
and  (18)  can  be  shown  to  be  consistent  with  eq  (18)* 

We  can  now  summarize  the  solution  proced-ure  for  the  principal  case. 

Note  first  that  by  the  Independence  Assumption^  the  constant  p . is 

^ J 

independent  of  i and  can  be  abbreviated  p.  = Of  where 

J J 

(19)  O'  = probability  that  a random  section  of  operator  lane  carries 
a parcel  to  be  sorted, 
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(*)  / 

TABLE  !♦  Values^  ^ ot 


=0.1 

Pj=0.2 

|Z= 

0 

1 

2 

3 

0 

1 

2 

3 

- 1 

•000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

2 

.100 

.001 

.000 

.000 

.200 

.012 

.001 

.000 

3 

.190 

.005 

.000 

.000 

*360 

.055 

.009 

.002 

4 

.271 

.014 

.001 

.000 

.488 

.154 

.054 

.020 

5 

.027 

.002 

.000 

.590 

.330 

.21+0 

.191 

6 

.410 

.049 

.006 

.001 

.672 

.589 

.703 

.788 

7 

.469 

.082 

.014 

.002 

• 738 

.868 

.995 

8 

.522 

.132 

.034 

.009 

.790 

.996 

9 

.570 

.208 

-091 

.041 

.832 

10 

.613 

.326 

.192 

.204 

.866 

11 

.652 

.502 

.710 

.735 

•893 

12 

.687 

♦735 

.985 

.914 

(*)  Omitted  values  are  i*00« 

Buffer  capacity  b(l_,j)=b  for  1 > 1« 
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TABLE  2 


• Values^  •'of  1-E.  , 
ij 


Pj=0 

.1 

Pj-0.2 

b= 

0 

1 

2 

3 

0 

1 

2 

3 

II 

H 

.000 

•000 

.000 

•000 

.000 

.000 

•OQO 

•QOO 

2 

.100 

.100 

.100 

*100 

.200 

.200 

•200 

.-200 

3 

.190 

.200 

.200 

.200 

.360 

•398 

•400 

•400 

ij- 

♦271 

.299 

.300 

.300 

.488 

t — 
00 
LT\ 

•597 

.600 

5 

.3^14 

• 399 

*400 

.i}00 

•590 

•756 

.786 

•796 

6 

.410 

.496 

.500 

.500 

.672 

.890 

•938 

•958 

7 

.469 

•591 

•599 

.600 

.738 

.971 

•997 

8 

.522 

.683 

•698 

.699 

*790 

.998 

9 

.970 

.770 

•794 

.799 

.832 

10 

•613 

•849 

.885 

.895 

•866 

11 

.652 

.916 

•969 

•973 

•893 

12 

.687 

.966 

.998 

— 

•914 

— 

_ — - 

•— 

FT 

Omitted 

values 

are  l.C 

)0* 

Buffer  capacity  b(l^j)=b  for  1 > 1* 
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(20)  f . = probability  that  a random  parcel  is  destined  for  the  j-th 
J 

cross  belto 

In  practice.  O'  might  be  interpreted  as  the  ratio  (if  < l)  of  operator 
coding  speed  in  parcels  / minute  to  operator  lane  speed  in  sections  / 
minute,  and  f^  as  the  relative  frequency  of  parcels  to  the  destinations 
of  the  j-th  cross  belt®  Th.e  initial  conditions  (for  l=l)  are 

(22)  = 1® 

For  the  recursion  step,  suppose  E,  „ and  P,  , are  known®  First 

1-1,  j 1-1,  J 

calculate  E.  . by  ttie  formula 
ij 


(23)  hj  = h-l,j  -ij  <•^-''1-1, 3^ <-''1-1 J 


b(l-l, j)+l, 


obtainable  from.  (’4)  and  (l8)®  fhen  compute  P.  . hj  (16)®  The  overflow 
rate  at  (i,j)  is  found  from  the  formula 

obtainable  from  (l),  (I7)  and  (l8)« 


This  procedure  was  used  to  obtain  the  numerical  results  presented 
in  Tables  1.  and  2®  The  tabulated  quantities  are 

V.  o/p.  = expected  fraction  of  parcels  on  the  i-th  operator  lane, 

destined  for  the  cross  belt,  which  overflow  at  (l,j). 


1 


= expected  fraction  of  the  cross  belt  sections  already 
occupied  when  they  reach  (i,j)« 


Although  the  table  refers  to  situations  in  which  the  buffer  capacity 
has  the  same  value  b(i,j)=b  for  all  1 > 1,  the  method  is  not  limited 
to  such  cases®  Additional  tables  are  in  preparation® 

In  discussing  these  numerical  results,  we  note  first  that  it  was 
obvious  from  the  start  that  adding  a.  buffer  (or  Increasing  buffer 
capacity)  at  an  "early"  Intersection  corresponding  to  a small  value  of  i. 


"I3  “ 


vMle  reducing  the  overflow  rate  at  that  intersection,  would  also 

tend  to  make  cross  belt  sections  leave  that  intersection  "loaded” 

more  often  than  before  and  thus  to  Increase  the  overflow  rate  at 

intersections  further  "downstream"*  It  is  not  too  surprising,  then, 

to  find  for  example  that  increasing  b from  1 to  2 diminishes  the 

overflow  rate  for  the  initial  lanes,  but  Increases  the  rate  for 

subsequent  lanes;  the  point  is  that  our  analysis  gives  precise 

information  as  to  which  lane  marks  the  beginning  of  this  "inversion" 

(1=11  for  p =Q*1,  1=6  for  p.=0*2)*  Also,  it  was  not  obvious  without 
J J 

specific  investigation  and  computation  that  for  fixed  b and  p^  the 

values  of  Instead  of  increasing  gradually  with  i,  would  Increase 

extremely  sharply  for  a certain  "critical  value"  of  1,  If  p^=Ool  and 

~~~~~~~~~~  J 

b=2  or  3 for  instance,  then  operator  lanes  past  the  tenth  will  contribute 

relatively  little  sorting  according  to  Table  1— -Table  2 shows  that 

cross  belt  sections  will  be  occupied  at  least  97^  of  the  time  when 

they  reach  such  lanes--- so  that  really  compelling  reasons  due  to  other 

factors  would  be  required  to  justify  a system  with  more  than  ten  lanes 

under  these  circumstances*  Similarly,  if  p,=0*2  and  b=2  or  3 then 

J 

the  addition  of  operator  lanes  past  the  fifth  seems  questionable*  These 
examples  should  Indicate  the  usefulness  of  dai^a  like  those  in  Tables 
1 and  2* 
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5o  SOLUTION  OF  THE  GENERAL  CASE;  APPLICATIONS 


We  return  now  to  the  general  analysis  interrupted  in  Section 

Recall  that  the  remaining  task  was  to  utilize  eqs  (6)  through  (9)  to 

determine  the  functions  in  terms  of  e(^^  and  For  this 

ijk  ij 

purpose  we  shall  prove^  by  induction  on  k^  that 


(25) 


ij,k+l  ij 


k 


= Z 

s=o 


[B 


(N+1) 

ijs 


bU9  + 

Ijs-'  ijk 


p'^Ai-eW] 

IJ  IJ 


for  o < k < b(i^j)o 

Note  that  under  the  conditions  (12)  of  the  principal  case^  the  first 
SUJH  on  the  right  vanishes  term  by  term  and  the  previous  eq  (15)  is 
regained;  this  shows  why  the  principal  case  admitted  so  especially 
elegant  a solutiono 

To  prove  eq  (25)^  we  first  observe  that  for  k=o  it  is  just  a 
rearrangement  of  eq  (6) » Next_,  for  o<k<b(i^j)  we  can  rewi*lte  eq  (7) 


as 


i,];k+l  ij 


[B 


7:^)  - b7?.]  + Bp)  Pp)[l-Ep)] 


ijk 


Ijk 


f (H)  (N)  (H) 

)®i.jk  Ilj  hj 


ijk  ^Ij  ^ ij 

, pP)  [1-eP>])  ; 

ij«k-l  -^1,1  '■  ij  ^ 


if  (25)  holds  with  k replaced  by  k-1^  then  the  term  [---]  in  the  last 
equation  is  equal  to 


k-1 

r 

s=o 


[BP^I) 


Bp)], 

IJS  IJS^ 


so  that  eq  (25)  holds  for  k as  well  as  k-l<>  Tnis  completes  the 
induction  proof® 

The  general  procedure  for  finding  the  functions  b(^^  can  be 

Ijk 

described  as  follows®  Eq  (25)  with  k=o^  yields  in  terms  of 
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and 


and 

ijo  ijo 


Therefore  it  also  yields  in  terms  of 

ijl  ijo 


(W+2)^  Next,  eq  (25)  with  k=l  yields  B^^^  in  terms  of  B^^^  . , 

ijo  ^ a \ j ^j2  ijl  ijl 

(W+l)  B^^^  ; by  the  last  two  sentences,  it  therefore  yields  an 

expression  for  B^^l  in  terms  of  B^^^  , and  ¥e  now 

lj2  IJQ  IJO  IJO 


continue  in  the  same  way;  for  o<k<b(i,j),  is  determined  as 

— — 1 Jk 

(N+k) 


. ^(W)  ^(N+l) 

an  expression  in  B.  . , B.  . , * « «,  B. . 

IJO'^  IJO  ^ ^ ijo 


After  this  is  done,  we 
substitute  the  resulting  expressions  into  eq  (9)  and  obtain  a difference 


equation  of  order  b(i,j)  for  the  function  B.  „ (t)  ^ , , , , . 

^ ^ ijo'  o Once  this  equation 

(with  appropriate  initial  conditions)  has  been  solved,  the  functions 

B ( "t- ) 

ijk^  ''  are  determined  by  the  expressions  already  found  for  them  in  terms 
of  B.  « (t) o 

ijo^ 

The  simplest  non-trivial  application  arises  for  b(i,j)=l,  loe.  »a 
buffer  with  room  for  only  one  parcel*  Here  eq  (25)  with  k=o  yields 

(26)  b'")  e(«  = - b!«]  + b(T  p(*”[1-e(T] 

Ijl  tLJ  IJ  IJO  IJO-'  IJO  ^IJ  IJ 


= - b(”)  [A”)  * pf”)  E?^]  ; 

IJO  IJO  "-^J  ■'^IJ  IJ 


if  this  is  solved  for  and  the  result  substituted  into  the  appropriate 


version 


IJO  ijl 


of  eq  (9)^  we  obtain  the  first-order  difference  equation 

. ......... 

ijo  ‘-•^ij  ij  ^j 


(27)  = q?^)  4-  b1^)  . ql^)(l-E^^'bl 


IJO 


IJ  IJ 


IJ 


- 16  - 


for  the  function  B^^^(t) 


(o)  _ 


If  the  time  t.  ' = t*,  was  sufficiently 
ij  iJ 

early  in  the  system's  operation,  we  may  take  = 1 as  initial 

conditlono  'fhe  solution  of  the '-difference  equation  is  then  found  to  he 

^ ^ 1.1  o ^,1  1,1  . id  ^ id 


-j  ij 
w-i 


s=o 


+ r 

s=o 


tj  hj 


W 

n 

r=s+l 


ij  ij 


where  the  last  sum  is  omitted  for  N=o;  this  formula  can  he  verified 
hy  induction  on  N« 

For  any  particular  "reasonable"  value  of  h(l_, d)  > ‘tke  difference 
equation  analogous  to  (27)  and  the  explicit  solution  analogous  to  (28) 
can  he  worked  out*  It  is  important  to  realize,  however,  that  for  many 
purposes  the  explicit  solution  is  unnecessary;  particular  values  of  it 
can  he  calcuJ,ated  from  the  initial  conditions  and  the  difference  equation, 
while  its  theoretical  properties  can  he  proved  hy  induction  with  the  aid 
of  the  difference  equation. 


As  a second  application,  let  us  see  how  much  the  analysis  of  the 
principal  case  is  complicated  if  (lO)  must  he  changed  to 


(29) 


k.  = 2 , 


l.e.  if  the  cross  helt  is  moving  twice  as  fast  as  the  operation  lanes. 
The  appropriate  modification  of  the  "steady  state  input"  condition 
(ll)  for  this  case  is 


(30) 


(2m) 

p:  . = P.  • j. 

ij  ±3  ^ 

(2m+l) 


Independent  of  m, 


since  we  automatically  have  p ^ = o when  (29)  holds.  Conditions 

(lOa)  are  retained,  and  the  Independence  Assumption  permits  us  to 

abbreviate  p. , to  p..  Under  these  hypotheses,  it  is  known  (see  the  end 
^ J «3 

of  this  section)  that  the  j-th.  cross  helt  will  gradually  settle  down 
so  that  the  analogs 


(31) 


E 


(W+2) 


,(N+2)  (N) 

'ijk  ^Ijk 
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of  (12)  are  satlsfiedo  (Eq  (12)  would  be  obtained  if  "N+2"  in  (31) 
were  replaced  by  "N+l",)  Thus  the  functions  E^j(t)  and  each 

take  on  at  most  two  different  values,  which  we  abbreviate 

E..  = Et,  = 

ij  ij  ^ iJ  iJ 


R _ R(2m)  + _ ^(2m+l) 

ijk  ^ijk  ^ ^ijk 


Eq  (25)  splits  into  two  equations; 

k 


(32) 

B.  . 

iJ^ 

I — 1 
+ 

H- 

C_1. 

II 

k 

(33) 

4- 

b;  . 

, ,1  eT.  = r 

We  now  define  the 

(3^) 

pij 

= 1 - 

(i/eT. 

■ iJ 

and  : 

prove 

that 

k 

[b|. 

(35) 

E 

- B.  . 

s=o 

IJS 

1 J s 

s=o 


IJS 


ijs^  ijk 


IJS  IJS' 


\ k + 
).  .)  B 
ij  o 


k 


s=o 


>..h-"  B.. 

ij-  IJS 


For  k=o,  eq  (35)  is  an  identity*  By  eqs  (33)  and  (3^)  f 
k+1  k 

r [bT.  - b.  , ] = r [bT.  - B.  . ] + bT.  , - B.  . , _ 

IJS  IJS  IJS  IJS  ij,k+l  ij,k+l 


k k 

= Z [bT,  - B.  . ] - (1/eT.)  E [B.  . - B.  , ] - B.  . 

“ IJS  ijs^  ' ij'  ^ IJS  ijs  ij,k+l 

S— O S— o 


= p.  . E [B.  . - B.  . ] - B.  . , T ; 

^ij  IJS  ijs-*  ij,k+l 


if  eq  (35)  is  true  for  k,  substitution  in  the  last  equation  shows  it 
to  be  true  for  k+1  as  well*  This  completes  the  Induction  proof* 
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The  case  h(i^j)=l  has  been  treated  earlier  (eqs  (27)  and  (28))^ 
while  the  situation  b(i_,j)=o  will  he  discussed  in  Section  6,  We 
therefore  assume  h(i_,j)=h  > 1,  and  proceed  to  show  that 


(36)  ®ij^k+l  ^^^ij  ^ij^  ^ij  " ^ j ^ 


B.  n.  . p.  (l-E.  ,)  for  o < k < b-2  ; 


ij 


this  is  an  analog  of  the  eg  (15)  obtained  for  the  principal  case, 
To  prove  (36)^  replace  k by  k+1  in  (32)  and  (35)  "to  obtain 

k+1 

B,.  . . ,p  q.  E = (p  ) B - r (p  ) B„  « 

ij^k+2  j ijo  ijs 


+ B.  . , p . [l-E.  . ] 
i^^k+1  J ^ 


p.  . [(p.  ,)^  B*.  -E  (p.  B.  . ] - B.  . , , 

^IJ  IJO  IJS^  ij^k+1 

0 — w 


k 


= p. . r [bT . - B.  . ] - B.  . , [g.  + p . E.  „ ] 

s=o  ijA+1  J J iJ 


^ij  ^®ij^k+l  \j  “ ®ijk  ” ®ij^k+l  ^ij 

yielding  (36)0  With  the  abbreviations 

(37)  A,  =._  p.  . (l-E,  ,)/g.E.  = ~p.  . V.  . > o 

ij  iJ  ^ J ij  iJ  - 

(38)  /i.  . = - [q.  + (p.=  - p.  . qJ  E }/g.  E , 

ij  J ij  J J iJ 

we  rewrite  eg  (36)  as 


(39)  \j  ®ljk  ^ ®ij;,k+l 
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This  second- order  linear  difference  equation  "with  constant  coefficients 
can  be  solved  by  standard  methods;  the  result  is 


hjk  = (hjl  - » hjo)  - y - 6) 

2 

where  y and  6 are  the  roots  of  the  equation  x = \.  . + fl.  . x » By  using 


b 


ij 


ij 


(40)  to  substitute  in  the  version  S B,  .,  = 1 of  eq_  (9)^  ve  obtain  the 


k=o 


ijk 


following  relation  between  B. . and  B. : 

ijo  ijl 

(in)  - (l-y)(l-6^’"^)} 

- C5(l-6)(l-y’"'"^)  - y (l-y)(l-6°*^)}  » (l-y)(l-6)(y-6). 

Another  relation  between  B.  . and  B.  will  be  derived  by  equating 
^ ijo  ijl  j-  a B 

two  expressions  for  To  find  the  first  expression,  we  set  k=o 

in  (32)  to  obtain 

(42)  bT.  = B.  q . E.  . + B,  . [q..  + p.  E.  . ] o 

^ ' ijo  ijl  ij  ijo 

To  find  the  second  expression,  we  set  k=b  in  ecq  (35)  and  use  the  facts 

= 1 

s=o  ^ s=o  ^ 


to  obtain 


, b , 

! \h  + « / \b-s  „ 

o = (n.  .)  B.  . - r (p.  .)  B.  . 

ijo  ijs 


or  equivalently 
b 

(43)  bT.  =1:  (p.  .)”®  B.  . ,* 

ijo  ijs 

use  of  eq  (40)  to  substitute  in  the  last  equation  yields 

<jo  = fhjl  [(«1J  - - (P,J  -y)((Aij)'*^^-6’="b'] 


^(Pij  -yKPij  -«)(Pij)’'  d-«) 
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The  relation  obtained  by  equating  the  right-hand  sides  of  eqs  (42) 
and  (44),  when  solved  simultaneously  with  (4l),  determine  and  B^.^. 

All  B^^^  are  then  determined  by  (40),  given  by  (42),  and  the 

remaining  bT ..  can  be  found  recursively  from  eq  (33) » This  completes 

ijh 

the  solution. 


We  conclude  this  section  by  sketching  how  (12)  follows  from  (lO) 
and  (lOa)  and  (ll),  and  how  (31)  follows  from  (29)  and  (lOa)  and  (30)* 
These  are  the  special  cases  k^=l,2  of  the  assertion  that  if  (lOa)  holds, 
and  if  the  Inputs  to  the  j-th  cross  belt  are  "steady  state"  in  the 


sense  that 


(w+K.)  (W) 

(45)  p.  . ,r  = p:  1 
ij  13 


for  all  W, 


then  the  cross  belt  will  gradually  settle  down  so  that  the  analogs 


(46) 

IJ  IJ  ijk  ijk 


of  (45)  are  satisfied. 

Number  the  "section  positions"  along  the  cross  belt  consecutively; 

the  first  position  is  at  (l,j),  the  (c.+l)-st  at  (2,j),  the  (2c  .*^-l)“St  at 

^ ^ ( c ) 

(3^j)j»  etc.  Define  the  state  of  the  system  at  time  t=t.+mk^t^  ^ 

J 3 3 

to  consist  of  the  number  of  parcels  in  each  buffer,  and  the 

status  (occupied  or  empty)  of  every  section  position,  at  time  t and  at 

each  of  the  preceding  moments.  It  can  then  be  shown  that  the 

probability  of  transition  from  any  particular  state  at  time  t to  any 

( c ) 

particular  state  at  time  t+k.t,  depends  only  on  the  states,  not  on  t. 
Thus  the  belt  constitutes  what  is  called  a Markov  chain,  and  the  assertion 
follows  by  appeal  to  the  theory  of  such  chains.  In  most  practical 
situations  the  "settling  down"  is  rapid  enough  that,  unless  especially 
interested  in  transient  effects,  one  can  confine  attention  to  the  problem 
as  characterized  by  the  equations  analogous  to  (46), 
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6,  THE  NO- BUFFER  CASE 


If  b(l,j) 
work  as 


then  the  fimctions  B.  (t) 

ijk 


are  known  without  further 


B.  . (t)  = 1,  B.  (t)  = o for  k > o« 
Ijo^  ' ^ ijk'  ' 


Therefore  It  is  relatively  easy  to  study  a cross  belt  with  no  buffers 
at  any  of  its  intersections,  or  merely  with  no  buffers  at  a set  of 
intersections  at  the  beginning  of  the  belt;  in  the  second  Instance  the 
investigation  is  confined  to  the  "buffer-free"  zone® 


For  this  situation,  eq  (5)  simplifies  to 


(47)  n 

J 


^,(47-  r = . t<7) 

r=a 


(i  > 1), 


and  so  eq  (l)  yields 

(48)  v7hp44i-  n 4..  (t!T-  S = . tUb] 

s=r 


(i  > 1) 


Suppose  we  Sspecialize  further  by  imposing  the  "homogeneity" 

(lOa),  Then  the  Independence  Assumption  shows  p. .(t)  and  q. .(t) 

1 J 1 J 

independent  of  i,  and  eq  (48)  becomes 

P^^^[1-  TT  q.  (t^.^^-  (i  - r+i)  c.t^.°^)]  (i  > l). 

ij  J ^^2  J J J 


conditions 
to  be 


o-  (n) 
Since  p.  . ■ 
ij 


o if  N is  not  divisible  by 


a non- zero  probability  of 


overflow  is  possible  only  in  the  case 

(49)  V..(t*  + mk.t^, ^^)  =p.(t*  + mk.t^, ^^)[l-  ff  Q..('t*  + [mk.  -(i-r+l)c.} 

ij  J 3 3 3 3 3 3 r=2  ^ ^ 3 3 


)] 


of  the  last  equation..  In  (49),  the  only  q. -values  not  necessarily  equal 

J 

to  unity  are  those  for  which  (i-r+l)c.  is  divisible  by  k . . If  for  example 

J J 

k.  and  c.  are  relatively  prime  (have  no  common  factor  except  unity),  then 
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the  only  non- unity  values  for  q.  in  (49)  are  those  for  which  i-r+1 

J 

is  a miiltiple  of  i.e. 


r=  i+l-k..  i+l-2k.,  i+l-3k..  etc.  > 2) 

J 3’  3’  ^ ^ 


The  number  of  such  values  is  < (i-l)/k.  > , where  <>  x > denotes  the 

J 

largest  integer  not  exceeding  x* 

Assume  now^  in  addition  to  (lOa),  that  k.  and  c.  are  relatively 

3 3 

prime  and  that  the  steady  state  condition  (45)  holds;  this  last  condition 
can  be  written 


p . (t*  + mk ,t ) = p, 
3 3 3 3 3 


independent  of  m_, 


and  (46)  permits  us  to  abbreviate 

fc) 

V. .(t*  + mk.t  \ M = V. independent  of  m. 

13  3 3 3 ij 

By  the  preceding  remarks^  the  overflow  probability  (49)  takes  the  form 

(50)  V = p (1-q (i“l)Aj  > ) (1  > 1). 

J J 


This  generalizes  a result  obtained  in  the  FMC  technical  memo  cited 

(see  "Step  pp.  J-8  of  the  memo)^  showing  a cross  belt  with  kj  > 1^ 

under  the  stated  hypotheses^  to  be  equivalent  to  several  shorter  belts 

each  with  k„=l.  Our  symbols  k.  and  c„  correspond  respectively  to  1 and  J 
J J J 

in  FMC ' s notation. 


If  we  specialize  still  further  to  the  case  k.=l 
and  let 


N=  number  of  operator  lanes 

(this  should  not  be  confused  with  our  earlier  use  of  the  symbol  W),  then 
it  follows  from  (50)  that  the  expected  rate  of  overflows  along  the  j-th 
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cross  belt  is 


N N N 

T V = r P,  (l-q,  ) = P,  (N-l-r  , 

x=z2  i=2  <J  cl  J i=2  ^ 


or  finally 


(51) 


N 

r 

i=2 


V.  . 

ij 


If  there  are  R cross  belts,  all  bufferless  and  with  k,=l,  and  if  all 

(c)  ^ 

c^  and  tj  are  equal,  then  the  expected  overflow  rate  for  the  entire 

system  is 


R N 

(52)  r T 


R 


R 


V. . = w r p.  - r (i-q^) 

j=l  i=2  J=1  ^ j=l  ^ 


Setting  p . = af 
J ' 
rate  becomes 


(see  eqs  (19)  and  (20)),  we  find  that  the  total  overflow 


(53)  0,  = cm  T f.  - T (i-(i-o:f  )^)* 

j=l  ^ J=1  ^ 

This  is  a generalization  of  the  corresponding  formula  on  p*  9 ol  "the 

FMC  memo,  which  can  be  obtained  by  setting  all  f ,=1/r  in  our  eq  (53)^  i«.eo 

J 

by  assuming  equal  distribution  of  parcels  among  cross  belts* 
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